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We present a complete precision analysis of the sfermion pair production process e^e fifj (/ = 
t, b, T, v-r) in the Minimal Supersymmetric Standard Model. Our results extend the previously calcu- 
lated weak corrections by including all one-loop corrections together with higher order QED correc- 
tions. We present the details of the analytical calculation and discuss the renormalization scheme. 
The numerical analysis shows the results for total cross-sections, forward-backward and left-right 
asymmetries. It is based on the SPSla' point from the SPA project. The complete corrections are 
^ about 10% and have to be taken into account in a high precision analysis. 

■ PACS numbers: 12.15.Lk, 12.60.Jv, 13.66.Hk, 14.80.Ly 
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I. INTRODUCTION 

^ '' The Minimal Supersymmetric Standard Model (MSSM) provides the most attractive extension of 

(~| , the Standard Model (SM). Among other particles it includes supersymmetric partners of the fermions. 

These scalar states /l, fa (sfermions) correspond to the two chirality states of each fermion /. The mass 
. eigenstates /i and /2 though are not identical with /i, fn and are rather a linear combination of them. 

/\ ' The mixing terms are proportional to the mass of the corresponding fermion. Hence the sfermions of the 

third generation play a special role. As a consequence, one eigenstate (/i) can be much lighter than the 
other one. 

The sfermions, especially the strongly interacting ones bi), are likely to be detected at the LHC or the 
Tevatron. Nevertheless, to extract the fundamental parameters one must have a significant accuracy only 
obtainable at a linear collider. From sfermion pair production in e"^e~ collisions the sfermion mixing angle 
can be extracted. This is one of the reasons why it has been extensively studied phenomenologically [1]. 
To match the expected precision of the linear collider, theory predictions must reach a similar accuracy. 
The effort to calculate higher order corrections to the sfermion production has begun by calculating 
the leading QCD, SUSY-QCD and Yukawa corrections [2-5]. It was further shown that taking only 
the leading terms of the one-loop corrections is not sufficient and so also the full weak corrections were 
presented in [6, 7]. 

It is the aim of this paper to extend the existing weak corrections by including the full 0{a) contributions 
in a similar manner as in the case of the selectrons and the smuons in [8] . In addition, we present the 
full analytical results and all the details of the calculation for both the weak corrections [6] and the 
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QED contributions. Moreover, we generalize the results to include also the effects of polarization of 
the electron and positron beams. Apart from cross-sections, we calculate other observables such as the 
forward-backward and the left-right asymmetries as well. 

Although we present the results in the form of cross-sections and asymmetries, we are well aware of 
the fact that the precise predictions have to be used for parameter extraction. As the definition of the 
parameters is no longer unique beyond the tree-level, there has been a recent proposal by the so-called 
SPA project (SUSY parameter analysis) which defines these parameters [9]. The SPA project also gives 
a firm base for calculating all sorts of observables (masses, decay widths, cross-sections etc.) and enables 
the development of tools for extracting the parameters. 

The fundamental SUSY parameters in the SPA project are defined using the DR (dimensional reduction) 
rcnormalization scheme at the scale Q = ITeV. Specifying the renormalization scheme serves only to 
define the parameters uniquely and does not restrict the use of other schemes in different calculations. 
In this paper, we use an on-shell renormalization scheme. To use the parameters from the SPA project 
we have to translate them into the on-shell renormalization scheme. The results for any observable using 
different schemes (with correctly translated input parameters) must agree up to contributions of higher 
order. 

The paper is organized as follows. In section II we give the formulae for the tree- level cross-section 
for polarized electron and positron beams. The calculation of the virtual corrections with a detailed 
discussion of the applied on-shell renormalization scheme arc outlined in section III. All explicit analytic 
formulas needed for the calculation arc given in the Appendices A, B, C and D. In section IV we work 

out the real radiative corrections where we include the Bremsstrahlung process a{e^e~ — > fifj^)- In 
section V we present the numerical analysis with some results of the corrections. Section VI summarizes 
our conclusions. 



II. TREE LEVEL 



The sfermion mixing is described by the diagonalization of the sfermion mass matrix given in the 
left-right basis (/^, fj^) into the mass basis (/i, /2), / = t, 6 or t [10], 




where is a 2 x 2 rotation matrix with rotation angle Bp which relates the mass eigenstates /, 
p < rrip) to the weak eigenstates /„, a = i, i?, b; 
= —^21 = smdf, and 



i = 1, 2, {nip < nip) to the weak eigenstates /„, a = i, i?, by /, = uf^fa, with = = cosOj and 



m?^ = M^^^^j -f (/|^-e/s^)cos2/3m| +m^, (2.2) 
m?^ = M'^^.Qj^^+efs'^cos2l3m^ + mj, (2.3) 
af = Af - fi {tan (3)-^''i\ (2.4) 



Mq, Mj^, Mjj, Mj^ and Afj, are soft SUSY breaking masses, Af is the trilinear scalar coupling parameter, 
H the higgsino mass parameter, tan /3 = ^ is the ratio of the vacuum expectation values of the two neutral 

Higgs doublet states , I^f^ denotes the third component of the weak isospin of the fermion /, e/ the electric 
charge in terms of the elementary charge e, and sw is the sine of the Weinberg angle 6w ■ 
The mass eigenvalues and the mixing angle are 

rn^flj = — {O<0f<7T), (2.6) 



and the mass of the sneutrino v-r is given by m?^ = 7\jf| + | cos 2/3. 

The tree-level cross-section of e^e~ — > for polarized electron and positron beams is given 

by 

a'-'' = \{l-P^){l + P+)a'r + \{l + P^){l-P+)at\ (2.7) 

where P_ , P+ G (—1,1) are the degrees of polarization of the electron and positron beams (e. g. P_ (P+ ) = 
—0.8 means 80% of electrons (positrons) left polarized and 20 % unpolarizcd). 

As we neglect the electron mass, we have only two terms contributing (out of 4 possible) where a^£'^'^ 

is the tree-level cross-sections for e^e£ — ^ Ji f^ (below referred to as the left part of the polarized cross- 
section) and cr^''° stands for eje^ fifj (analogously referred to as the right part of the cross-section). 
They have the form 

^t%e+e- ^ fj^ = tf {tZ^^ + T2% + TE%) , (2.8) 



where 



e 



^^'^ " 32(s-m|)2 ^^■^^> 



.zz _ 9zio-ij) ^2 
and k{x, y, z) = ■\/(x — y — z)^ — 4yz. 

Here we use i^L,/?, and Cl,/? as the left- and right-handed couplings of the electron to the photon and 
Z-boson, respectively, 

Kl = Kr = 1, Cl = -\ + sI^, Cr = sI,. (2.12) 



The matrix elements a- come from the coupling of Zfifj^ 



, _ , 4(/f cos2 0;-.^e/) -2/fsin2(?; 



2/|^ sin 26 f A{Pf^ cos^ 6l^~ - s^Cf) 



Apart from the tree-level cross-section we can calculate other observables such as the left-right asymmetry 
and the forward-backward asymmetry. They are defined by 

A,n - , = ^^-^ , (2.14) 

with 

^ ^ d^j^ ' (^) dcos^, I (^) rfcos.. (2.15) 

There is no lowest order (tree- level) contribution to the Aj?B-asymmetry as the angle distribution is 
symmetric. 



III. VIRTUAL CORRECTIONS 



For a precision analysis of the sfcrmion production one has to include also higher order corrections. 
The calculation of the higher order corrections is performed analytically in the DR scheme, adopting the 
^ = 1 'tHooft-Feynman gauge. All necessary ingredients of the analytical calculation are given in the 
Appendices. Furthermore, we neglect the electron mass wherever possible (toc = 0). For the numerical 
evaluation of the loop integrals we use the packages LoopTools and FF [11]. At the end the whole analytic 
result was checked with the result obtained using the computer algebra tools FeynArts and FormCalc 

[12]- ..... 

The virtual corrections receive contributions from vertex, self-energy and box diagrams depicted generi- 
cally in Fig. 1 and explicitly in Figs. 2 and 3. All these contributions are summarized in the renormalized 
cross-section o-''°". 

The one-loop (renormalized) cross-section 0'''°" for polarized beams is expressed analogously to Eq. (2.7), 

^""(e+e- - /J,) = \{l P-){1 + P+) + J(l + ^-)(1 - P+) <^T > (3-1) 
where the left/right renormalized cross-sections are defined as 

_.ron ^trco , a QCD , \ _EW fq n\ 

with the symbol A denoting UV-finite quantities. 

The SUSY-QCD corrections (Acr^'-'^) have already been calculated for the unpolarized case in [3, 4]. As 
the gluon part of /S.a'^'^^ is proportional to the tree-level cross-section, the polarized cross-sections arc 
easily obtained using cr*^'^^ instead of cr*'''^*^. The gluino part of A(t*3CD treated analogously to /S.a^^ 
(see Fig. 1). 
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FIG. 1: UV-finite parts of the electro- weak corrections 

We have already presented the unpolarized results for the electro- weak corrections (Act^^) in [6]. In this 
paper, we give the result for polarized beams and also all formulas needed for the calculation. 
The electroweak corrections can be split further into four UV-finite parts given in Fig. 1, 

Aaf% ^ AaY% + AaY% + A<?? + A<« , (3-3) 

where Ac^^'J^ and Acr^'^ stand for the left/right part of the renormalized electron and sfcrmion vertex, 
A(j^'^°j^ and Acr^°^ for the left/right part of renormalized propagators and box contribution. 
The renormalized electron vertex has the form 

= f^ ^'^^'^r^^ ((AT^r/).,. + (A^r/),,, + (ATJI),^ J , (3.4) 

where 



{AT^,')^^^ = {^eL.RKL,R), (3.5) 



(^'ryz)L..R = - Zis-ml) (^^^'^^^.^ + ^^lmKl.r) , (3.6) 



(AT^i)^^ = ^M4^(Aai,^Ci,H). (3.7) 

16 (S — TO^j 



AcL^fi and Aa^^/j consist of 3 parts, 



^eL,n = Se[% + Se[-l + Se^L^H. (3-8) 
AaL.R - -Ja^^^.A^ -^ai'A + ^"l fl ■ (3-9) 



'^^i'fl' correspond to the vertex corrections in Fig. 2, (^e^^]j, (^a^"']j are the wave-function corrections, 

and Se^^^j^, Sa^f'^j^ correspond to the countcrtcrms. 
The rcnormahzed sfermion vertex has a similar form, 

Nr 



^-1% = ^ ^ AtV ^ „ + AtV , „ + AT-/ . J , (3.10) 



where 



(Ar,vX^ . SiSi^^Al,^, (3,11) 



r^-^^ = - 4,(;i^|) Kl.mCl,r {{Aef),,ai^ + %(Aa/),,) , (3.12) 



at;/ _ = ^S^^Ci,. (3.13) 



16 (s — m|)2 

(Ae/)y and (Aa/)ij can also be split into vertex corrections (see Fig. 2), wave-fmrction corrections and 
counterterms, 



[Aefh = {5ef)f^ + {Scfp + {5ef)'t;> , (3.14) 
{Aaf\, = {5aft^ + {5a}p + {5ait;>. (3.15) 



The diagrams contributing to the vertex corrections are shown in Fig. 2 and the explicit form of the 
corrections are given in Appendix A. The wave-function corrections and the counterterms to both ver- 
tices are listed in detail in sections IIIAl and III A 2. The {5ef)ij and {5af)ij corresponding to gluino 
corrections can be found in [4]. 

The correction AaY°ii which comes from inserting the self-energies of the 7 and Z-boson, in the propa- 
gator, see Fig. 3, can be expressed as 

AapP = ^ 




swcw 



UtDl.r + {Th)L,B) + ^ {iT,%)L,R + {T^z)l.,r 



where T2}j, T^^^, Tf_| are defined in Eqs. (2.9-2.11) and 

^^-rii^'R - 4s(s-m|) 2 



[J^Z'yJL.R - -~-r-r„ Z:TTn^L,R^ (3.17j 




FIG. 2: Vertex diagrams contributing to ^e^"^ , 5al\i , {Sef)]^ and ((5a/)'j . 



^^^^^^'^ = 16(s-l|)2 2^^'^^'^^'«' ^^-^^^ 

{T^.y)L,R = -^Cl,rKl,r , (3.19) 

^^^^^^'^ ~ ~4s(.-m|) 2^^-«- (^■^'^^ 

The nyy(s) in Eq. (3.16) are the transverse parts of the renormaUzed self-energies of the vector bosons 

7 and Z. The unrcnormahzed self-energies are given in Appendix C2 and the renormalization is done 
following [13]. 

The box corrections arc obtained by adding up the diagrams shown in Fig. 3 and are given by 

Aa^°S = ^ T^V^ / sin^ ^ d^ , (3.21) 



where 



4 4 7rs"^ 



^tfl - I - (4^)2 2^'^'^ + (4^)2 4(r- m|) 2^^'^ ' ^^'^ ' ^^'^^^ 

The Bl,/? are the form-factors defined in Appendix B, where one can find the analytic expressions as 
well. 





A. Renormalization scheme 



In order to make the result finite we have to introduce the wave-function renormahzation constants 
and counterterms. We fix them following the on-shell renormalization scheme. The parameters already 
occurring in the Standard Model (SM) are renormalized according to [13]. We assume the CKM matrix 
to be diagonal and so have no fiavour mixing among the SM fermions at one-loop level. 



1. Wave-function renormalization 

The wave-function corrections are due to a shift from unrenormalized (bare) fields to the renormalized 
(physical) ones. For the fields relevant here we have 



Z^) \ \5Zz^ 1 + \5Zzz )\Zt, 
The form of the corrections for the left vertex is 



(3.23) 



Jei^-jj = {5ZLM. + ¥Z^i)KL.R-Yf5Zz^CL,R. (3-24) 

5a^i^\ = {5Zlm. + ¥Zzz)Cl.r-\—5Z.,zKl.r- (3.25) 

^ gz 

The wave-function corrections for the right vertex are 

(•^e/)!]"^ = \ef{5Z,j+5Zji) + ]-ef5Z^^5ij + —^ a{~- (5Zz^ , (3.26) 

= ^'^[SZk^al^+6ZkJa{A+2swcwef6iJSZ-^z + ^a{J5Zzz, (3.27) 



where sw = sinOw and cw = cos 6*1^ • 

The wave-function renormahzation constants are determined by imposing the on-shell renormaUzation 
conditions as in [14, 15] such that the on-shell masses are the real parts of the poles of the propagator 
and the fields are properly normalized, 



2 5R4(m|) 



SZ,, = 



SZ, 

SZzz = -^^zz{mz)-, ^^iz = - 
SZl 



2 9 



2 5Rn^z(m|) 



SZz-y = 



-g?ri^-y(o), 

2 3finz^(0) 



= 



5Zf 



-ULiml) - mlinLiml) + URimD) 
nieitlsLiml) +tlsR{ml)) , 

nniml) - mlitlniml) + tlLiml)) 
meitlsRiml) + nsL{ml)) 



2m, 



-{UsLimD^nsRiml)) 



(3.28) 
(3.29) 

(3.30) 



1 



2mf 



(nsRimD-UsUml)) 



(3.31) 

where ri(m^) = [■^Il{k'^)'\ We use the self-energies given in Appendix C and in [16] where we 

adopted the conventions from. 

A remark should be made at this point. Wc include the wave- functions renormahzation constants of 
the vector bosons although they arc not external particles. By introducing them into the wave-function 
renormahzation of the vertices, we have additional checks that can be made. First of all, the renormal- 
ization constants of the vector bosons must drop out in the final result. Secondly, the vertex corrections 
and the propagators can be both made UV-finite separately. 



2. Counterterms 



The counterterms come from the shift from the bare to the physical parameters in the lagrangian. It 
includes the shifting of e, mw, 'fnz^ 0^ defined by 



mw + Smw , = mz + Smz , 



The counterterm contributions for both vertices are 



Se 



Sa 



L,R 
L,B. 



6e 

- Kl,r , 
e 



Se 
e 



5 m 



w 



Smz \ 1 — 2sw f Smw Sm^ 



mw mz 



J \ mw mz 



(3.32) 

(3.33) 
(3.34) 



('5e/)!? 
('5«/)!f 



— e/ 

e 



/ Smw Smz 
\ mw mz 



2 f ^"^W 

ScfCw I 



Sm2 



\ mw mz 



(3.35) 
(3.36) 



where the contributions containing S9^ were intentionally left out and will be discussed below. 



3. Renormalization of electric charge 



The standard on-shcU input value for the electric charge is the one in the Thomson limit a = e^/(47r) 
1/137.036. This corresponds to a countcrterm 



Se 
e 



2cw 



5Z 



(3.37) 



In this way of fixing the electric charge has a significant theoretical uncertainty coming from the 
light quarks which we circumvent by using as input parameter for a the MS value at the Z-pole, 
a = C({mz)\^ = e^/(47r). The counterterm then is given by [16-18] 



Se 
e 



(47r)2 6 



4 5: A^/ e?.(A + log J) + t iV/ e?.(A + log ^) 
/ ^ ^ S ™=i 

+ 4EfA + logJ 



fc=i 



A + log^ 



21(A + log^ 



(3.38) 



with Xf = mz y rrif < mz and xt = mt- n[, is the color factor, = 1,3 for (s)leptons and (s)quarks, 
respectively. A denotes the UV divergence factor, A = 2/e — 7 + log47r. 



4- Renormalization of mw and mz 

The masses of the Z-boson and the M^-boson are fixed as the physical (pole) masses, i. e. 
6ml = 5ini^(m|), 

where 

Smz 1 Sm^ 
mz 2 ' 

The formulas for the vector boson self-energies ^^^y(m^) and Hl^{mD are given in Appendix C and 
in [16]. The counterterms for the intermediate boson masses are used to determine the Weinberg angle 
fixing according to [19]. 



(5m^ = 5Rn^,^(m^), (3.39) 



Sm^^lSm^^ (3.40) 
mw 2 TO^ 



5. Renormalization of 6? 



The countcrterm of the sfermion mixing angle. 60j, is fixed such that it cancels the anti-hermitian part 
of the sfermion wave-function corrections [5, 20], 

% - \ iSZ,2-SZ,,) = \ ^ ^ {nUml)+uUrn}j) ■ (3.41) 

^ fi h> 

Including the terms proportional to 56j in Eq. (3.36) is equivalent to symmetrizing the off-diagonal 
sfermion wave-function corrections in Eq. (3.28) as [15, 21] 

5Rn{2(m2 )-5Rn{i(m2 ) 

5Z^2 = 5Z21 = 2 —■ (3-42) 

m % — m% 

Jl J2 

This fixing of the counterterm for the mixing angle is analogous to the renormalization of the CKM 
matrix in [22] and similarly has to be made gauge- independent. It was shown in [23] that this can be 



avoided or, equivalently, the result in the ^ = 1 gauge can be regarded as the gauge-independent one. 



IV. REAL PHOTON CORRECTIONS 



Similarly to the QCD case where the cross-section was IR-divergcnt due to massless gluons [2-4, 24], 

the one-loop cross-section (7''°"(e+e^ fifj) is IR-divergent owing to the diagrams with photon exchange 
where the photon mass is zero. This is remedied by introducing a small mass A and including also the 
Bremsstrahlung process i. e. (T(e+e^ — * fifjj), see Fig. 4. Summing these two contributions yields an 
IR-finite result for the physical value A = 0, 

a™"(e+e- ^ fj^) = a-"(e+e- ^ fj^) + a(e+e- ^ fj^-f) . (4.1) 

To calculate the radiative cross-section a{e^e~ — *■ fifjj) we use the phase-space splicing method [25] 
which splits the bremsstrahlung phase-space into 3 regions. The corresponding 3 parts are 

a{e+e- fjjj) cr"°f*(A, A^;) + o-hard^^^^ a''°^\AE, A9) . (4.2) 

In our calculation, we used a soft-photon approximation ((t^°^') to reproduce the divergence pattern 
correctly. However, this approximation introduces a cut AE on the energy of the radiated photon. The 
dependence on the cut AE drops out if we include the full 2^3 process {a^'^'^'^). In order to get simpler 
expressions for a^'^'^'^ we neglect the electron mass but then a collinear divergence occurs when the photon 
is radiated in the direction of the electron and positron beams. This collinear divergence can be regulated 
by introducing yet another approximation (cr'^"'') for the above mentioned phase-space region. Another 
cut A6 is hereby introduced. After summing the 3 contributions the result must be independent of both 
the cuts and has to cancel the IR-divergcnce of the one-loop cross-section. This is the ultimate test we 
have made at the end of the calculation. 




A. Soft-photon approximation 



The soft-photon approximation supposes that the 4-momentum of the photon is small compared to 
other momenta (for details see e. g.[13]). Using this assumption the differential cross-section (f^j^^fj is 
proportional to the tree-level differential cross-section. The full cross-section for polarized beams is 



a^°'\\,AE) = i(l-P_)(l + P+X* + i(l + P_)(l-P+)a|?ft, (4.3) 



where 



soft 



^"^'^^ dn^f(^] 6.dn. (4.4) 



/soft J \ dn y^^^^ 



The factor Ss is defined as 

Os = -- 



^ {^pj + ^pI ^ 2^PiP2 + 6/(4^ + Ikl - + 2ef{Ip,ki + Ip2k2 - Ipik2 - Ip2ki)j ,(4-5) 

where the integrals lat are defined in [13] and were worked out e. g. in [26]. The expHcit formula for Ss 
can be found in Appendix D. 

B. Hard and coUinear photon radiation 

The cross-section for the full brcmsstrahlung process e+(p2) e~(pi) fi{ki) fj{k2) 7(^3) is given by 
ah-d(A£;,A0) = Ys8(hy J \M\'dk°dkUvdcose, (4.6) 

where the cuts AE and A9 appear in the integration bounds of dk^ and dcos6'. The angle 77 is defined 
as in [12]. The explicit form of the squared matrix element is given in Appendix D and the integral is 
evaluated numerically using the routines from the CUBA library [27] . 

As we have neglected the electron mass in the calculation of a^^^'^, we have to take an another approach 
in the coUinear region of the phase-space. We follow the approach of [25, 28] and get for the coUinear 
cross-section the following expression , 

a™'>(Ai?,A0) = l[(l-P_)(l-P+X' + (l-P_)(l + P+)a™A' 

+ (1 + P-)(l - P+) ati^ + (1 + P_)(l + P+) a'j^j^] , (4.7) 

where all polarization states {o-f^^, cr'j^j^, c^x', cT/fj?') appear. This is due to the radiation of an additional 
photon and the fact that the electron is massive {aLR stands for the cross-section with left-handed 
electrons and right-handed positrons etc.). The single polarization states are given by 

<L,RR = ^ '^^ ^ - + ((1 ~ x)s)] , (4.8) 



<^TIhl = ^ I dx2 



'L,R 



{{l-x)s), (4.9) 



with 



(mf + m ? )^ 

Xrmn = 2AE/^, Xraa. = 1 . (4.10) 

s 

After including all the above-mentioned contributions we arrive at a cut-independent result. 

C. Higher order corrections 

Substantial correction from the coUinear photon radiation is due to the smallness of the electron mass 
compared to a typical energy scale in the process. This effect is such that to reach the collider precision 
one has to include also the leading higher order corrections (i. e. beyond 0{a)). Owing to the mass- 
factorization theorem, one can factorize the corrections in the leading-log (LL) approximation as 



y"da"™ + y"d(7LL = ^ dxij^ da;2rLL(xi,Q')rLL(a;2,Q2) Jda'''^%x,p,,X2P2) 



(4.11) 



where Xi, X2 are the momentum fractions of the electron and the positron carried after the radiation of 
the photon(s). 



The (a;, Q^) is the leadmg-log structure function up to 0{a'^), given in ref.[29], 



exp(-l/37i; + f/3) /3 
r(l + f) 2 

-^(l + ^) + Y^(-2(l + x)log(l-x)-^ + ^(l+x)logx-|-^) 



-a + x) 



32-^ + 4log(l-.)+2log^(l-x)-- 



loga-log(l — x) 



1,2 1twi ^^ ll + a;2 / 3, 1,2 1, , .1 n\ 

-log -Li2(l - I + -y— -(^- -log.T+ — log X- -logxlog(l-x)J 

-i(l-a-)(log(l-x) + i) +^(5-3x)logx- , (4.12) 



with the gamma function F, the Euler constant ~ 0.577216, and /3 = ^(log^ — 1). For the free 

scale we take the typical energy of the process s. The soft-photon contributions were summed up to 
all orders in the perturbation scries. 

The structure function (4.12) contains not only the higher orders beginning with 0{a^) but also parts of 
terms 0{a) already included elsewhere. To avoid double counting we subtract these terms as in [28]. 



V. NUMERICAL ANALYSIS 



In contrast to [6] , we do not attempt to make a scan over a large area of MSSM parameter space but 
rather consider only one benchmark point in the numerical analysis. It is the SPSla' point we use as 
input which is defined in the SPA project [9]. The point is chosen such that it satisfies all the precision 
data and both the bounds for the masses of the SUSY particles and the bounds from cosmology. 
The input parameters for the SPSla' point are defined in the DR scheme at the scale Q = ITeV. As 
we use the on-shell renormalization scheme, we have to transform the SPSla' input parameters V into 
on-shell parameters 7^*-*^. This transformation is simply performed by subtracting the corresponding 
countertcrms i. e. V'^'^ = "PiQ) ^ SV{Q) and the results for the relevant parameters are listed in Table 1. 
All other parameters do not enter in the calculation at tree-level and so the differences when using 
the on-shell or the DR value are of a higher order. A further remark is necessary here. One of the 
parameters not entering the tree-level directly is the infamous Ai, parameter. Fortunately, the parameter 
set taken here causes all the on-shell input parameters to be insensitive to the problems of the Ab 
on-shell definition. 



V 


DR 


OS 




M 


OS 




402.87 


399.94 






368.6 


tan/3 


10 


10.31 






583.1 


Mq., 


471.26 


507.23 






499.9 


Mo, 


384.59 


410.11 




mr 

02 


543.7 


Mn, 


501.35 


538.92 






107.5 


Me, 


109.87 


111.58 






195.4 


Ml, 


179.49 


181.78 






170.7 



TABLE I: Parameters of the SPSla' scenario in the DR and on-shell scheme and particle masses. 



The Figs. 5-8 show the total cross-sections for the pair production of the sfermions of the third 
generation. In general, we show the complete corrections and the tree-level where the tree-level is 
defined according to the SPA project. According to the SPA project, all masses arc taken on-shell and 



all parameters in the couplings are given in the DR scheme. The virtue of using this tree-level definition 
is that not only the total corrected cross-sections are directly comparable to other calculations using the 
SPA conventions but one can also compare the relative corrections. 
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FIG. 5: Left; The total cross-sections for stop pair production (all channels) at tree-level {dashed} and with 
complete corrections {solid}. Right: The total cross-sections for polarized beams: stands for P_ — —0.8, P+ = 
0.6 and a+ for P- = 0.8, P+ = -0.6. 




FIG. 6: Left: The total cross-sections for sbottom pair production (all channels) at tree-level {dashed} and with 
complete corrections {solid}. Right: The total cross-sections for polarized beams stands for P_ = —0.8, P+ = 
0.6 and a+ for P_ = 0.8, P+ = -0.6. 

For each sfermion type we show an unpolarized case (left) and a case where the beams are polarized 
(right). We take two sets of polarizations, either P_ = —0.8 and P+ = 0.6 or P_ = 0.8 and P+ = —0.6. 
The difference to the earlier calculations [6, 7] are the QED corrections which give a negative contribution 
near the threshold due to the known soft-photon behaviour. The QED corrections are substantial (as can 
be checked when comparing the results of this paper with those of [6]) and cannot to be neglected. The 
plots on the right-hand side of Fig. 5-8 show the effect of beams polarization on the radiative corrections. 
Polarization and its effects are best seen in other observables which we discuss in the following. 
The Figs. 9-12 show the left-right and forward-backward asymmetries for different final states as 
defined in Eq. (2.14). Owing to the fact that at tree level there is only a s-channcl contribution the 
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FIG. 8: Left: The total cross-sections for tau-sneutrino pair production at tree-level {dashed} and with complete 
corrections {solid}. Right: The total cross-sections for polarized beams a~ stands for P- = —0.8, P+ = 0.6 and 
0-+ for P- = 0.8, P+ = -0.6. 



■\/s-dependence drops out in the left-right asymmetry, making it to a good approximation constant. The 
•\/s-dependence is then a result of the one-loop corrections. Notice that the corrections are substantial 
especially in the tit2, '?i'?2: as well as in the Vr'^T channel, where there is only a Z exchange at 

tree-level. 

As we have already mentioned, there is no tree-level contribution to the forward-backward asymmetry 
and thus the asymmetry is loop-induced. In the calculation of the forward-backward asymmetry at 
one-loop one has to define the forward direction, in particular for the contributions coming from the 
photon radiation . We define it by ap = cr(cos0-,- > 0) where 9^,~ is the angle between the 
incoming electron and the outgoing sfermion with negative isospin. As an additional feature, we also 
show the forward-backward asymmetry for polarized beams where the polarizations are P-_ = —0.8 and 
P+ = 0.6. In general, one sees that the asymmetries receive sizeable corrections and thus justify the 
higher-order calculation. 
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FIG. 9: Left: The left-right asymmetry for stop pair production at tree-level {dashed} and with complete cor- 
rections {solid}. Right: The left-right asymmetry for sbottom pair production at tree-level {dashed} and with 
complete corrections {solid}. 
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FIG. 10: Left: The left-right asymmetry for stau pair production at tree-level {dashed} and with complete 
corrections {solid}. Right: The left-right asymmetry for tau-sneutrino pair production at tree- level {dashed} and 
with complete corrections {solid}. 



VI. CONCLUSION 



We have calculated the full 0{a) corrections to stop, sbottom, stau and tau-sneutrino production in 
the MSSM. We have presented the details of our analytical calculation which was also checked by the 
computer algebra tools FeynArts and FormCalc [12]. The results extend our previous calculations [4-6] 
by including also QED contributions and the real photon radiation. We have also used the structure 
function approach [28] to include some higher order effects. Moreover, the whole calculation was 
extended to the case of polarized e^-beams. 

In the numerical analysis, we have studied only one specific scenario based on the SPSla' benchmark 
point defined in the SPA project. We have transformed the input parameters into the on-shcll 
renormalization scheme which we have used throughout the paper. The numerical results show the total 
cross-sections and asymmetries with the effect of the 0{a) corrections. These are found to be sizeable 





FIG. 12: Left: The forward-backward asymmetry for sbottom pair production (no tree-level contribution) for 
unpolarized beams. Right: The forward-backward asymmetry for sbottom pair production with polarized beams 
P- = -0.8, P+ = 0.6. 



(in some cases up to 15% and larger), and in particular the forward-backward asymmetry is only due to 
higher order corrections. 
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APPENDIX A: VERTEX CORRECTIONS 



Here we give the explicit form of the clcctrowcak contributions to the vertex corrections which are 
depicted in Fig. 2. For SUSY-QCD contributions we refer to [4]. Ah couphngs used in this paper can be 
found in [16]. 

The vertex corrections iSef)[^^ and (<5a/)|j'' (or <5e^"^^ and (Ja^"^^) originate from the diagrams in 
Fig. 2 with 7, Z-boson exchange, respectively. 

The vertex corrections to the right vertex (5e/)|j'' and {Saf)[^'^ are defined on the amplitude level from 
the corresponding diagrams. The general form of the amplitude is given by 



= al^2 , HP2)riY^ Pl + yX Pr) u{p,) [A,, ih - fe)^ + By (fci + h)^] . (Al) 



For Xy and Y^j^ with V" = 7, Z we have 

r7 



X^=e, YIj^ = Kl,r, Xz^-gz, Y^^j, = Cl,r. (A2) 

We use the form-factor Aij (the other form-factor vanishes in |A^P) to define the vertex correction iSef)[j^ 
and (5a/)|j' as 



The explicit formulas for (<5e/)^J'' and {5af)'f^ are given below in the Appendices Al and A2. 

The vertex corrections to the left vertex, (5e^^^ and include corrections to the two chiral parts of 

the electron photon/Z-boson vertex. The generic form is 



iXv)i;j 

(47r)2 s — my 

+ClPl{pi + P2Y + CrPr{pi + p2Y]u{pi){ki ~ k2)f, , (A4) 



= 773^2 ^(P2) [^L7^ Pl + Pr + BlPl{pi - P2T + BrPr{p^ - ps) 



where {Xv)ij {V = 7, Z) stands for 



gz f 

{X^),j = -e efSij , {Xzh = — ^af^ ■ (A5) 

From this generic structure only the form-factors Al^r survive. We define the vertex corrections as 



The contributions to the vertex corrections ^ef^^i ^'^l^r ^^"6 given in the Appendices A 3 and A 4 



1. Corrections to "yfifj vertex 
The vertex correction {6ef)[^'^ is composed of contributions from different classes of diagrams as follows, 

(^^/)? = (^^/)?'*^ + iSejf'' + ('^e/)!;-^' + (<5e,)i;-^+^' + (Se,)^^^^ + (Sej)^^^^ . (A7) 

In the following we use the standard two- and three-point functions Bi and Ci from [30] in the conventions 
of [13]. We introduce the following standard set of arguments Ci = Ciirn^- , s,my, Mq, M'^ , M2) to be 

used in the generic functions S. 



The first contribution coming from the exchange of one or two gauginos is 

with the generic vertex function 

S,f^{Mo,A'h,M2;g§,gl^,g^,g^,g^,g^) = M,M2{g^g{gk + g^g^g^){Co + Ci + C2) 

+MoMi(go^gf52« + g^g^g^) (Co + Ci + C2) + A'/iM2(5o'^<7f Ss"- + 9o9i92) (Ci + C2) 

+ {g^gh2 +9o9f9k){Bo{s,MlMl)+Ml{2Co + Ci+C2)+m)C^+m)C2) . (A9) 

The corrections due to graphs with 3 scalar particles in the loop are given by 

iSef)[f^ = -T^'^lf Y.G{LG^MC, + C,+C2){ml,s,ml,m),^,rn\^,ml^) 

^ ' k,m=l ^ k k 

^ 2 _ _ 

- ^/ ^LkGLjk (Co + C1+C2) {ml , s, m| , m^o , m^^^^ , mj^ ) 

-7i^^/' E G£G{£,(Co + Ci + C2)(m|,.,m|,m^+,m2-^,m2^). (AlO) 

^ ' k,m=l 

The graphs with one vector particle (7, W^) and two sfermions in the loop yield 

I 2 _ _ 2 ^ 
+74;;ti ye/' i?,fii?Ji ^ 5,^^^(mvi/,TO/^,mj, ) , (All) 

where SY/H- . .) is a short form for 

S^ff{Mo, Ml, M2) = 4Coo + (m2 - 2s + 2m| + 2m|) (Co + Ci + C2) 

+ (2m| + 2m2.^ ~ s) (Ci + C2 + Cn + 2Ci2 + C22) + (™| - ) (Ci - C2 + Cn - C22) ■ 

(A12) 

From the diagrams with one W^+-boson, one Goldstone boson and a sfermion we obtain 

1 2 

cAp \i".ssv+svs} _ I gmw sr^ {pf rff' 4- rff) nf' 

V / V 771 — 1 

X (Co - Ci - C2)(m|,s,m^^ , TOy:, ,m^,TO^) , (A13) 



and with the generic scalar-vector-vector vertex function 

(Mo, A'/i, M2) = 2Bo {s, Ml M|) - 2Coo + + SM^) Co - s{Ci + C2) 

~i(m| - m|.) (Ci -C2 + Cu - C22) - i(2m| + 2m| ~ s) {Cn + 2Ci2 + C22) (A14) 
the correction due to the exchange of two W^+-bosons and one sferniion reads 

i^^ft"^^ - 7A1 E {RL)"s^{m^-,^,m^,m^) . (A15) 

^ ^ ra—l 

The contributions from one sfermion and one vector particle (7, Z", W~^) can be expressed as 

1 



1 2 „ ^ 



rn— 1 



1 2 _ _ ^ 

-7^ ^ E (^mi)'(2So + Si) (m^^, , m^,, m}^) + ^ j) (A16) 

with y/ = 2(/^^ — e/). The symbol (i <-> j) denotes the previous term with the indices i and j 
interchanged. 

2. Corrections to Z^fifj vertex 

The corrections to the Z'^f^fj vertex have the same components as in Eq. (A7). Using the same 
abbreviations for the generic vertex functions as in the previous section we get for the single contributions: 



4 



k,l=l 

2 



(An) 



+ 7Iv[2 E ^^^^ {"^xl ' ^i'^L^L ajfe' 



.2 

■^^5f^(m^+,m/sm/,;C^\c{\4,;t;/,,fc/,) , (A17) 
ype sferm: 

sfermions (down-squarks and sleptons) 

9-2 4 

iSaf)^u''''^ = ^ Y^R,,_2{P-a)G{^,Gi^,{Co + C, + C2){ml,s,rnl,mj^,m],o,mlo) 

^ ' k,m=l 1=3 

" E E Ri^k-2{P - a) Gl^^Giji {C0 + G1+C2) , s, m| , m^^ , m^o , m^o) 



with Oj^^^ = Oj^^^ for up-type sfermions (up-squarks and sneutrinos) and O^i^^ = O^^^ for down-type 



4 



(47r) 



4 



/f cos 2ew Gi£,GfCk (Co + Ci + C2) (™| , m| , m}^ , m^^+ , m^+) 



fc,m— 1 

2 



^LnGirnkGijk (^0 + Cj + C2) (m| , s, m^. , m^o , , m^.^) 



(47r) 

/c.m.n— 1 



2 



^ E ^LG{£kGili{Co + Ci+C2){ml,s,ml,ml^,mj^,mj^), (A18) 



(47r) 

A;,m,n — 1 



^ ^ ' m,n— 1 

^ E i?£ii?,{;4^^~(miv,m^^, ,m;,), (A19) 



(47r)2 2 

m,n— 1 



/■- -.(v.SSV+SVS) 

('5a/ )V 



(47r)^ -^-^ V -^"^ "''jl""2m2y ""mi V^u ^ ) V fi'> ^ f^^ f!^ 

^ rn— 1 

+ 7^^ 2 E i?fe2 (/3 - «) 4mGlnjk (^o - (^i - C2) , s, m| , to|^ , to|, m^o) 

^ ' fc,m=l 

J 2 _ _ 

"'"TiTrP '^9zmz ^ i?*;^ (/3 - a) z{^,jG{^^ (Cq - Ci - C2) (to| , s, m| , to|^ , m^o , m|) , (A20) 

^ ' k,m=l ' 

i^-ff^^ - yAlS^^^^f 4i4 E (i?,(:i)'5,f^(™/;^,-v^,-w) , (A21) 

^ ' m=l 

4 2 ^ ^ 



-, ^ ^ 2 _ 

--2 yI rIr(, E (^:™i)' (2^0 + Si) (m| , m^^, m|^) + {z ^ j) . (A22) 



(4^) 



3. Corrections to 'ye'^e vertex 

In the following wc list the analytic formulas of the vertex corrections to the electron-positron-photon 
vertex. We only give the right-handed coefficients of the generic vertex functions, ■* as the coefficients 

can be obtained by exchanging the indices R and L, i. e. ^ = \r L). 
In the remaining vertex corrections we use the standard set of arguments for the whole class of C-functions 
C = C{mlsM,MlMlM^). 

The vertex correction Se'j^^-p^ is split into the following classes: 

-^e^ = ^eff « + ,5e(,-^« + Set^' + ^^t^' (A23) 



The contribution from one sfermion and two gauginos in the loop is given by 

2 



k=l 

where we have used the generic vertex function 

-{g^h'^'h^" - 9^9k9kM^)Co - {g^gM'' ~ gUkh\^)^^C, + {g^g^h^^ - g^gih?^)m,C2 

(A25) 

and the abbreviations (no sum over i) hf'^ = g{me + gfMi for i = 1,2 and (j, k) = L, R. 
The corrections due to the exchange of one gaugino and two sfermions are 

Se^n'^"^ = tAiE E {bl,yC,o{rnl,s,rnl,mlo,rnl^,mlJ , (A26) 

' k=lm=l 

Using the generic vertex function for one vector particle and two fermions in the loop, 

J'X^^{Mo,A'h,M2;gl^,gl^,g^,gtgk,g2) = ~ 2[gl^ g^ (2C00 - Bo(s, Af^M^) - ml{Ci + C2) 

-M^Co + (s - 2ml) {Co + Ci + C2) + + g^g^g^M.Ah C'o - g^g^g^m^ {Co + Ci + C2)] , 

(A27) 

with r = in the DR renormalization scheme we get the corrections stemming from one vector boson 
(7, Z") and two electrons given by 



5e 



{v,VFF) 
L,R 



j^e" T^il (a, me, me; 1, 1, 1, 1, 1, l) 
+ j^glT^!'i[mz,m,,m,-A,l,C'k,Cl,C'j,,Cl) . (A28) 
For the graphs with one electron-neutrino and two M^-bosons we obtain 

^^tr"^ = ^^^i;^(0,m^,m^;l,0,l,0,l), (A29) 
where we have used the function 

Ti^{Mo, Afi, M2; go, gt 91,92. 92) = 5o [5iV(2So(s, Af^, Af|) - r + 4Coo + 2A/2Co 

+ (5m2 - 2s) (Ci + C2)) + i{gig^ + g?gk)m,MoCo + igigk^UCi + ^2)! . (A30) 



4. Corrections to Z^e'^e vertex 

In the following we list the single contributions to the electron-positron-Z*^ vertex. The generic vertex 
functions used in this section can be looked up in Appendix A3. 

^ i 2 

^'^L.^R^^^ = ~ (/^^y E E -^LJi f "^e„. , m^o , m^ ; O^./^ , 0^/" , Qmfc , Knk . Kil 7 ^ml ) 
^ ' k,l=lm=l 



(47r) 



^ 2 



k,l=l 
4 2 



^ ' k=lm,n=l 

1^-2 

"•"74^2 5Z ('^ifc) '^oo('7ie:S,mg,m|o,m?^,m?J , (A32) 

^ k=l ^ 

and Ja^"'*^"-^"^ = Sa^^;'^^ ^\bl^ af., fc^^ ^ 



1 



Jai'^'r''^ = TT^e^ (a, me, me; C^,C£, 1,1, 1,1) 



9% ^Lji {j^z, Trie, rrie; C^, C£, Cj^, C£, C^, C£^ 

^ ^ •?iT(™>^' 0' 0; 0' 1' 0' 1' 0' 1) ' (A33) 



(47r)2 
1 

Sair^ = --i^^^i:^(0,m^,mH.; 1,0, 1,0,1). (A34) 



(47r)2 2 

APPENDIX B: BOX CONTRIBUTIONS 

In this section we give the exphcit form of the radiative corrections which stem from box diagrams 
with two different topologies. The matrix element is parameterized as 



<Pi), (Bl) 



where the form-factors Al r do not contribute to the squared matrix element. The form-factors Bl r 
depend on the Mandelstam variables t and u which are defined as 



t = -^(m^ +m^f - s) + I; K{s,m']: ^m"] ) cos d , (B2) 
ji J 3 

The single contributions to the form-factors 



u = 2ml + "^i + ^A. ~ s - t. (B3) 



Blm = Bl\ + Bl' + i?f,| + B^^ + BU + BU (B4) 



correspond to the diagrams with two vector bosons, where the particles in the loop arc indicated by a 
superscript, and similarly B^ ^ and Bf^ ^ denote the contributions from charginos and ncutralinos, 
respectively. 

1. Vector bosons in the loop 

In the case of two vector bosons, we use the generic functions 

{Mo, Ml, M2, M3) = 4Co + C2 + iM^Do + M^Ds +4[{mj^~u)Di - sDi + (t-m|)Z?2 + t-Da] 

(B5) 



for a diagram with vector bosons V in the loop and 

(Mo, Ml, M2, M3) = - {iCS + Gl + ^MlDl + ^{u-m^~^^ )Df + 4{u~mj}D^ + {M^ + 'iu)D^^ 

(B6) 

for the corresponding crossed diagram. In case of a M^-boson there is no crossed diagram and we use 

or depending on the charge of the final state particle. 

The scalar three-point and four-point functions used above have a standard set of arguments defined as 



a = a{s,mj^,ml,MlMlMi) 
Cf = a{mj^,ml,s,MlMi,Ml) 



Df = D,{mls,m] ,u,mlm] ,M§,MlMlMl) . 



The contributions from 2 photons, one photon and one Z-boson, 2 Z-bosons and 2 H^-bosons are 



n77 



e^e)5^, { B 



2 



(me, 0, 0,771 , 



nZZ 



efdz^L.R ^ Sirnzl^j B^ [nie, 0, Ulz , m jf^) -f B^^(777e, 7772, 0,777^-^ ) 
m— 1 
2 

-e^efg%CL.R ^ 4ni^^3 B^ {rne,mz,Q,m^J) + B^"" {me,0,mz,m^^) 

m—l 

2 



Br 



(B8) 

(B9) 
(BIO) 
(Bll) 



where B^^^^ denotes B^^ for up-type sfermions and B^ for down-type sfermions in the final state. 

2. Scalars and fermions in the loop 

In analogy to the case with vector bosons in the loop we define the following generic function for box 
diagrams with fermions and sfermions in the loop 



[Mo, Ml, M2, Afa; 9(^, <?f , 9^,92 , 92, 9z, <?3 ) = 
9o9i 9293 (Co{s,m)^,m\,MlMlMl) + MlDo - {Ml^t)D^ 
-PhM2 gl^g^g^giiDo + D-s) - M1M3 glig^g^g^D^ - M2M3 g^g^g^g^D^ 
and for the crossed counterpart 

Bl^^Mo, M,,M2,M^; g^, g^,g^, gt g^ , g^,gi, g^) = 

9o9i 9293 {Coiml , vfi\,s, MI,MI,MI) + MID^ ~ {Mi-u)D^) 



with B 



FF{x) 
R. 



B 



FF{x) 



{L ^ R). 



As in the case of two ly-bosons in the loop, for the graphs with charginos we use either 



(B12) 



M1M2 giigt92 93 iD^o + D3) ~ M^Nh gi^gtg^gSDI - M2M3 g^g'^g^g^D^ , (B13) 



(B14) 



for up-type sfermions or 

2 

^IT = E ^rHKe,™;f^«;^^m/.;fcf^,Zr^,4,ziz|,fc/,,/rF,40 (B15) 
for down-type sfermions in the final state. 

The contributions from 2 neutralinos in the loop have the following explicit form: 

^442 

fe=l 1=1 m=l 
4 4 2 

fc=l 1=1 m=l 

APPENDIX C: SELF-ENERGIES 

Here we give the explicit form of the self-energies needed for the computation of some wave-function 
renormalization constants and various counterterms. We omit the sfermion self-energies already given in 
[16]. All fermion. sfermion and vector self-energy diagrams are shown in Figs. 3 and 13. 

1. Fermion self-energies 

In our notation, the fermion self-energy is defined as 




M = iu{k)n{k)u{k) 



with 

u{k) = ^PLn^(fc')+ /^Pfl,n«(fc2) + n'5^(fc2)Pi + n^«(fc2)p^. (ci) 

Below we list the contributions to the left- and right-handed parts H^^^ and H^l.sr ^youi the single 
diagrams. The form- factor 11 is defined as a sum of the contributions coming from the diagrams in 
Fig. 13. 

n = W"" + n^-^t + U^i" + + W + + W"^* . (C2) 

We give the full formulas for the electron self-energy without neglecting the electron mass (although it is 
being neglected in the actual calculation). 

Note that for quarks and leptons (contrary to charginos), the left- and right-handed scalar parts of n(fc) 
are equal, i. e. U'^^ik) = U^^{k). 
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FIG. 13: Fermion and sfermion self-energies 



2. Vector self-energies 

Here we give the explicit form of the general gauge boson self-energies (the transverse parts only) which 
are then applied to the cases of the photon and the Z-boson (and their mixing). The corresponding 
couplings are given in a table after each generic formula. We do not list the contributions to the 
counterterms of the Z- and M^-bosons as they can be found in [16]. 
The self-energy of a vector boson is defined as follows 




The transverse part of the self-energy consists of the following parts: 



+{ulyr +(,nlyy-'^— (CIO) 

For the contributions with a fcrmion loop we define the foUowing generic function: 

Avv{k\mi,m2;gt , 92,92) = 2 [ - {g^g^ + g^g^) {k^B^ + mjBo + Ao{ml) - 2Boo) 

+ {9i92 +9i92)mim2Bo\ (k^mlml) (Cll) 
Using the generic function we can write the 2 fcrmion, 2 neutrahno and the 2 chargino contributions as 
{U^y)"{k') = j^Y.^C^^v{k',mf,mf;g^,g?,g^,g^), (C12) 
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TABLE 11: Couplings for the 2 fermion, 2 neutralino and the 2 chargino contributions to 77, 7Z and ZZ self- 
energies. 

The next set of contributions are the ones with 2 scalar particles in the loop. For this set we introduce 



Bvv{k^,ml,ml;gi,g2) = -4.9132^00(^^7711,7712) 
and get for the sfcrmion, neutral and charged Higgs in the loop the following forms: 

1 



n 



T \ff f i„2 



VV ) 



(fc^) 



(47r)2 



X! ^vv{k'^,m^f,rn^fjgi,g2) 



/ m.n—l 
2 4 



(4.)2 ^ 

{^iv)"^"hk') = 7iE«vv(^^-^.,r7l^. 

^ ' fc=l 



; .91 > 92 ) 



(C15) 

(C16) 
(C17) 
(C18) 



The next class are the self-energies with a single scalar particle in the loop for which we use the generic 
form 



Cvv{m^]gi) = giAoini^) 



(C19) 
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TABLE III: Couplings for the 2 sfermion, 2 neutral and 2 charged Higgs contributions to 77, 7Z and ZZ self- 
energies. 



The diagrams with 1 sfermion, 1 neutral or charged boson can be written as 
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TABLE IV: Couplings for the 1 sfermion, 1 neutral and charged Higgs contributions to 77, 7Z and ZZ self- 
energies. 



The diagrams with a vector and a scalar particle in the loop use the simple generic form 

'Dvv{k'^,ml,ml;gi,g2) = gig2 Bo{k'^ ,ml,ml) (C23) 

and give 
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The remaining 3 contributions comprising of 2 M^-bosons, 2 FP ghosts and a single VF-boson in the loop 
have the following explicit forms: 
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TABLE V: Couplings for the Z-boson-neutral Higgs and the W^-boson-charged Higgs contributions to 77, 
and ZZ self-energies. 
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TABLE VL Couplings for the 2 H^-bosons, a single VK-boson and 2 FP ghosts contributions to 77, ^Z and 
self-energies. 



APPENDIX D: BREMSSTRAHLUNG INTEGRALS 
1. Soft photon integral 

Using the kinematics of the process we get for the explicit form 
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where 5, t and u arc the Mandelstam variables, d = Cj zb ^i: with q and defined as 



s + m| — . 



(Dl) 



(D2) 



2. Hard photon integrals 

The squared matrix element for the hard photon radiation can be spUt into 3 parts, 



(D3) 



where jA^-'^p stands for the part of the amphtude where the photon is radiated off the particle indicated. 
The squared matrix part corresponding to the photon being radiated from the electron or positron has 
the form 



IM^ = -(1 - + P+) \M^\l + -(1 + - P+) \M% . 



The chiral L, R parts are 
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with Srcd = [pi +P2-hY- 

The functions TZi and 7?.2 contain only scalar products of the external momenta and are defined as 
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+2to| (fcg.ps) + 2m| (fc3.p2) - 4 (fci./cs) (fc3.p2) 



(D9) 



[2 (fci.fcs) (fci.pi) + 2 (fci./cg) (fci.ps) - 4 (fci.pi) (fci.ps) - 2 (fci.pi) {k2.h) 



(pi./c3)(p2-fc3) 
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The radiation off the sfcrmion can be written as 
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(Dll) 



where 
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The tensor T^^ is defined as 
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The interference term of the squared hard photon amphtude is 

2n{M^M''^) = i(l-P_)(l + P+)23fJ(7W^A1°t)i + i(i + p_)(i-P+)25R(7W^A1°t)j 
The chiral L, R parts are 
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with Sj-cd = (Pi +P2 - fcs)^- 

The functions TZ^ and 7?.4 are given by 
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